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ABSTRACT 

We have used the SubmiUimeter Array to image a flux hmited sample of seven submilhmcter galaxies, 
selected by the AzTEC camera on the JCMT at 1.1 mm, in the COSMOS field at 890 /xm with 
~ 2" resolution. All of the sources - two radio-bright and five radio-dim - are detected as single 
point-sources at high significance (> 6cr), with positions accurate to 0.2" that enable counterpart 
identification at other wavelengths observed with similarly high angular resolution. All seven have 
IRAC counterparts, but only two have secure counterparts in deep HST/ACS imaging. As compared 
to the two radio-bright sources in the sample, and those in previous studies, the five radio-dim sources 
in the sample (1) have systematically higher submillimeter-to-radio flux ratios, (2) have lower IRAC 
3.6-8.0 /im fluxes, and (3) are not detected at 24/im. These properties, combined with size constraints 
at 890 /zm {9 ^ 1.2"), suggest that the radio-dim submillimeter galaxies represent a population of 
very dusty starbursts, with physical scales similar to local ultraluminous infrared galaxies, and an 
average redshift higher than radio-bright sources. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
starburst - galaxies: submillimeter - galaxies: formation 



1. INTRODUCTION 

Early studies of the far-infrared (FIR) cosmic back- 
ground indicated that up to half of the cosmic energy 
density is gener ated by dusty star b ursts and activ e 
galactic nuclei (jFixsen et al.l Il998t iPei et al.l I1999D . 
One of the most exciting developements of the past 
decade has been the resolution of a significant fraction 
of this background into discrete so urces. Deep , wide 
blank-field s urveys at 850 /zm (ISmail et al.l 119971: 
Barger et all 119981: iHughes et al.l 19981: lEales et all 
19991 120001: iCowie et al.l [200^ IScott et all [20F 



User Bolometric Array (SCUBA; iHolland et all Il999( l 
on the James Clerk Maxwell Telescope ( JCMT), and 
later s urveys at millimet er wav el engths (jCreve et al.l 



Webb et a l.' 2003: 'S erjeant et al.1 '2003': W ang et ai1l20& 
Coppin et al. ,2006il with the Submillimeter Common 



^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 
^ jyounger@cfa.harvard.edu 

^ Astronomy Department, University of Massachusetts, 
Amherst, MA, 01003 

Joint ALMA Office, El Golf 40, Las Condes, Santiago 7550108, 
Chile 

^ National Astronomical Observatory of Japan (NAOJ), 2-21-1 
Osawa, Mitaka, Tokyo 181-8588 

® Institute of Astronomy, University of Tokyo, 2-21—1, Osawa, 
Mitaka, Tokyo 181-0015 

Nobeyama Radio Observatory, Minamimaki, Minamisaku, 
Nagano 384-1805 

* Instituto Nacional de Astroffsica, Optica y Electronica 
(INAOE), Tonantzintla, Peubla, Mexico 

^ Department of Physics, McGill University, Rutherford Physics 
Building, 3600 rue University, Montreal, Canada, H3A 2T8 

Max-Planck— Institut fiir Astronomic, Koningstuhl— 17, D- 
69117, Heidelberg, Germany 

Astronomy and Space Sciences Department, Sejong Univer- 
sity, 98 Kwangjin-gu, Kunja-dong, Seoul, 143-747, Korea 

Astronomy Department, Smith College, Clark Science Cen- 
ter, Northampton, MA 01060 

Fellow of the International Max Planck Research School for 
Astronomy and Cosmic Physics 



20041: iDan ncrbau er et all 12004 ': Carilh et al." '200^ 
Schlaerth et al. 20051 : iLaurent et al ...2005,; .Bertoldi et al 
20071) ■ revealed that this background was dominated 
by luminous (LIRG) and ultraluminous (ULIRG) in- 
frared galaxies at high redshift Multi-wavelength 
follow-up studies of these submillimeter galaxies (SMGs) 
showed that they are massive, young objects seen during 
their formation epoch, with very high specific star 
formation rates that may account for up to ~ 50% 
of the cosmic st ar formation at z > 1 (see review by 
iBlain et al.ll200l ). 

Progress towards a thorough understanding of the 
physical processes driving SMGs has been hampered 
by two factors: their faintness at optical wavelengths, 
and the relatively poor ( ^ 10") angular resolution of 
the current generation of submillimeter cameras. The 
first significant breakthrough came with deep radio sur- 
veys, which found a correlation be tween subm i llime- 
ter and radio contin u um emission ftvison e t al.l 119981: 
Chapman et al.l 120011: llvison et al.l ^002; Du nlop et all 



20041: llvison et al.ll2007f) . This locahzed SMGs to a few 
tenth's of an arcsec ond, and allowed the first s pectro- 
scopic observations (jChapman et al.l l2003l l2005f l which 
showed that SMGs lie at high redshift 2 z <^'i with a 
median of z ^ 2.5. However, only a fr action of these red- 
shifts have been confirmed via CO (Greve et al.' '200| 
Tacconi et a l. 2006) or mi d-infrared (Lutz ct al. 200 
Menendez-Delmestre et al.l 120071: lyahante et al.l 1200 



spectroscopy. Furthermore, the rapid dimming of the 
radi o continuum with redshift (/ (1 + z)^^'^^°'\ a = 
0.8: lCondonlll99l ) means existing radio-confirmed SMC 
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samples, whi ch represent ^ 3/4 of the overaU SMG pop- 
ulation (e.g.. lIvison et al.ll2002D . are relatively insensitive 
to systems at z ^ 3, and thus are biased. Recent studies 
have suggested that near to mid-i nfrared imaging u sing 
the Infrared Array Camera (IRAC: [F^zio et alJ[2004[ ). in 
combination with 24 /zm observat ions using the M ulti- 
band Imaging Photometer fMIPSr lRieke et al.|[2004l ). on 
board the Spitzer Space Telescope may offer an alterna- 
tive to radio identification. However, t his technique relies 
on either purely statistical argument s (iPope e t al. 20061 
or broad band infrared color criteria (|Ashbv e t al. 2006 ) , 
and also may be subject to biases. 

Reliable counterpart identification represents poten- 
tially the most challenging obstacle to a more complete 
understanding of SMGs. Previous interferom etric obser- 
vations at millimet er (Downes ct al. 1999; Fraver et al\ 
2000tlDannerbauer e t al. 2002; Downes & Solomon 2003; 



Genzel et all 20031: iKneib et al.ii2005t iGreve e t al. 2005; 
Tacconi et al.ll2006l) and submillimeter ( lono e t al. 2006) 



wavelengths have identified unambiguous counterparts 
for increasing numbers of radio-detected SMGs, and 
have confirmed the radio-submillimeter association. 
However, to date there has been no reliable high- 
resolution foUowup of a uniformly selected sample in- 
cluding radio-undetected SMGs, and the true nature of 
these sources therefore remains elusive. 

In this work, we present high-resolution 890 fim inter- 
ferometric ima ging by the Submillimeter Array (SMA: 
iHo et all l2004( ) of a flux-limited sample of sources se- 
lect ed at 1.1 mm by the AzTEC Camera (jWilson et al.l 
[200^ on th e JCMT, in a surv ey of a section of the COS- 
MOS field (jScott et al.ll2007D . The SMA has confirmed 
all seven of the AzTEC targets at arcsecond resolution, 
with positions accurate to 0.2". In §[2] we describe our 
observations, and in §[3] we address sources of uncertainty 
in the derived positions. In §|4]we describe each of the 
sources, and in §[5] we discuss some potential interpreta- 
tions of the data. All magnitudes are given in the AB 
system (|Okelll974D . 

2. OBSERVATIONS AND DATA REDUCTION 



The COSMOS field f Scoville et all '2006) benefits from 
an extraordinary wealth of deep, multi wavelength cov- 
erage from the X-ray to the radio. In this work, 
we utilize i band imagin g with the Advanced Cam- 
era for Surveys (ACS: Fo rd et all [1998) on board the 
Hubble Space Telescope to a depth of 27.1 magnitudes 
(jKoekemoer et al.ll2007l ) , a variety of gr ound-based opti- 
cal and near-infrared im aging data (see iTaniguchi et al.l 
[20061 iCaoaket al.l[200l . IRAC and MIPS imaging at 
3.6, 4.5, 5.8, 8.0, and 24 ^m to 5cr depths of ~0 .9, IJ , 
11.3, 14.6, and 71 ^Jy respectively ( Sanders et al.ll2007[ ). 
and 1.4 GHz radio continuum imaging to a mean rms 
depth of ~ 10.5 //Jy/beam w ith the Very Large Array 
(VLA: iSchinnerer etal][200^ . 

The AzTEC/COSMOS survey covers 0.15 deg^ of the 
COSMOS fi eld at 1.1 mm wit h an rms noise level of 1.3 
mJy/beam (jScott et al.ll200l . The AzTEC/COSMOS 
catalog includes 44 sources with S/N > 3.5a and 10 ro- 
bust sources with S/N > 5a. For our SMA observations 
we chose the seven highest significance sources, which 
effectively yielded a flux-limited sample of millimeter se- 
lected SMGs. 

Five of these sources (AzTECl-4 and AzTEC6) have 



either weak (J20cm < 60 /LtJy) or no radio sources within 
the AzTEC beam (18" FWHM), which we designate 
radio-dim. The remaining two sources (AzTEC5 and 
AzTEC7), with strong radio sources {F20cm = 161 & 
196 /iJy) within the AzTEC beam are designated radio- 
bright. This convention was chosen to address the inher- 
ent ambiguity of the radio detected versus undetected 
design ation often used in the literature; as iPope et al] 
(|2006( l observed, SMGs without a radio counterpart 
likely do not represent a distinct population, but rather 
lie just below the detection threshold for a given survey. 

The SMA observations were performed in the com- 
pact array configuration (beam size ~ 2") at 345 GHz 
(full bandwidth 2 GHz) from January through March 
2007. The weather was excellent, with typical rms noise 
levels of 1.0-1.5 mJy per track with ~ 6 hours of on- 
source integration. T he data were calibra ted using the 
MIR software package (jScoville et al.|[l993l ). modified for 
the SMA. Complex gain calibration was performed using 
the calibrator sources J1058-I-015 (~ 3 Jy, ~ 15° away 
from targets) and J0854-t-201 (~ 1 Jy, ~ 24° away from 
targets). Passband calibration was done using available 
strong calibrator sources, primarily 3C273 and Callisto. 
The absolute fiux scale was set using observations of Cal- 
listo and is estimated to be accurate to better than 20%. 
Positions and fiuxes of the COSMOS sources were de- 
rived from the c a.librated visibiliti es using the miriad 
software package (jSault et al.lll995( ). 

3. ASTROMETRIC UNCERTAINTIES FROM 
INTERFEROMETRIC IMAGING 

Precise astrometry is one of the most valuable contri- 
butions of interferometric observations to the study of 
SMGs, and accurate characterization of the positional 
uncertainty is crucial. There are two factors to con- 
sider when estimating astrometric accuracy of SMA ob- 
servations of SMGs: (1) statistical errors due to noise in 
fitting a point-source to the calibrated visibilities, and 
(2) systematic errors due t o uncertainties in th e interfer- 
ometer baselines (see e.g.. lDownes et allll999[ ). In gen- 
eral, the 1-D statistical uncertainty in position scales 
as ~ 0.59/{S/N), where (S/N) is the sign al-to-noise 
of th e fit and is the FWHM of the beam (Reid et all 
Il988f l. This expectation is borne out in the MIRIAD fit- 
ting routines, which for the ~ 2" FWHM SMA beam 
and (S/N) « 10 yield typical uncertainties of ~ 0.1" in 
a and S. Systematic uncertainties, or errors related to 
uncertainties in the baselines, scale as ~ A{As/X)Rd, 
where As is the baseline error, R is the distance from 
the calibrator with known position in radians, and A is 
a constant of order unity that is sensitive to the details 
of the array-source geometry. For the SMA compact ar- 
ray configuration, the baseline parameters are typically 
measured to better than 0.1 millimeters rms. To obtain 
an empirical upper limit on the systematic position error 
induced by baseline errors, we use one of the calibrators, 
J1058-h015, to calibrate the other, J0854-h201 (35 de- 
grees away, more than twice the distance of J1058+015 to 
the COSMOS field), and examine the resulting offsets.^"' 
This procedure yields a typical systematic total angular 
offset of J0854-I-201 of ^ 0.2" from its known position. 

-^■^ The calibrator sources are sufficiently strong that statistical 
errors are insignificant. 



3 



TABLE 1 

AsTROMETRY OF SMA/AzTEC Sources 





Name 


a{a) 




AzTEC Offset 


IRAC Offset" 


AzTECl 


AzTEC J095942.86+022938.2 


0.11" 


0.20" 


3.3" 


0.3" 


AzTEC2 


AzTEC J100008.05+022612.2 


0.13" 


0.23" 


0.3" 


... t> 


AzTECS 


AzTEC J100020.70+023520.5 


0.19" 


0.31" 


1.6" 


0.9" 


AzTEC4 


AzTEC J095931. 72+023044.0 


0.15" 


0.24" 


3.5" 


0.5" 


AzTEC5 


AzTEC J100019.75+023204.4 


0.16" 


0.11" 


1.7" 


0.8" 


AzTECe 


AzTEC J100006.50+023837.7 


0.19" 


0.28" 


2.8" 


0.7" 


AzTECT 


AzTEC J100018.06+024830.5 


0.24" 


0.29" 


1.5" 


0.5" 



^ Relative to IRAC Channel 1 source, which has an astrometric uncertainty of ~ 0.2" and 
an angular resolution of 1.6". 

^ IRAC counterpart is confused with a bright foreground object. 



We combine this conservative estimate of the systematic 
uncertainties with the measured statistical error, in the 
uncertainties hstcd in Table [TJ 

4. NOTES ON INDIVIDUAL OBJECTS 

Astrometry and photometry for all the targets are 
included in Table [1] and [5] respectively, and postage 
stamps are shown in Figure [TJ Here we comment on the 
individual objects. 

AzTEC J095942.9+022938.2 (AzTECl) ~ AzTECl 
is the brightest submillimeter source in our sample, 
and is detected at 14a significance by the SMA. A 
weak radio source (F20cm = 48 ± 14 fiJy) is coincident 
with the SMA position. There are also IRAC 3.6, 4.5, 
and 8.0 /xm sources coincident with the SMA position, 
but no significant MIPS 24 /im emission. There is a 
compact B-band dropout offset from the SMA position 
by 0.1" with i = 25.25 ± 0.79 mag in the ACS mosaic 
which we believe is the optical counterpart. The 
B-band dropout nature of this so urce suggests that 
it lies at 3.5 <, z ^ 4.5 (see e.g., ISteidel et aD 119991 : 
iGiavalisco et al.1 l2004f ). which is consistent with the 
lack of strong radio or 24 fun emission expected from a 
starburst galaxy at that redshift. 

AzTEC J100008. 0+022612.2 (AzTEC2) ~ AzTEC2 
is detected at 12a significance by the SMA. There is 
a weak radio source (F20cm = 52 ± 14 /iJy) coincident 
with the SMA position, but it has no candidate optical 
counterpart. ACS imaging reveals that it is offset by 
~ 3" from a foreground galaxy that is very bright in 
IRAC and MIPS. Even a careful subtraction, using the 
ACS data convolved with IRAC and MIPS point spread 
functions to remove the foreground object, does not 
reveal a potential counterpart; AzTEC2 is either not 
detected in the near and mid infrared, or is severely 
confused with a foreground galaxy that is not associated 
with the submillimeter continuum emission. 

AzTEC J100020.7+023520.5 (AzTEC3) - AzTEC3 is 
detected at 6a by the SMA. There are IRAC detections 
at 3.6 and 4.5 fim coincident with the SMA position, 
but no likely MIPS 24 /im or radio counterparts. There 
is an optical detection with i — 25.91 ± 1.07 in the ACS 
mosaic within 0.3" of the SMA position. The lack of 
24 /im emission suggests that the source is at z ^ 3, 
or is at z ~ 1.5 and has a deep rest-frame 9.7 /im 



silicate absorption feature. From the obs erved 1.1 mm 
flux, and assuming a grey-body (see lYun fc Carillil 
[2003 ) with Td « 40 - 50K and » 1.5 - 2 and the 
local FIR-radio correlation of ICondonI ()1992[ ). we would 
expect strong 20cm counterparts with F20cm 100 /.tJy 
at z ^ 1.5. Therefore if this optical source is at 
lower redshift, it either has lower radio emission than 
would be expected from the local FIR-radio relation, 
or is a chance alignment and not the correct counterpart. 

AzTEC J095931. 7+023044.0 (AzTEC4) - AzTEC4 is 
detected at 7a by the SMA. It has no candidate optical, 
radio, or 24 /im counterparts, but is detected by IRAC 
at 3.6, 4.5, 5.8, and S.Qfim. 

AzTEC J100019. 8+023204.2 (AzTEC5) - AzTEC5 is 
detected at 8a by the SMA. It has a radio counterpart 
coincident with the SMA position, and is detected by 
IRAC at 3.6, 4.5, 5.8, and 8.0/im, and by MIPS at 
24/im. There are two significant radio sources within 
the AzTEC beam, with fluxes of F20cm = 161 ± 35 and 
81 ± 12 /iJy. One of these radio sources is singled out 
as the counterpart by the high angular resolution SMA 
imaging. Furthermore, the correct 24 /im counterpart is 
the weaker of the two within the AzTEC beam. There 
are, however, no associated optical sources in the ACS 
i band image. The IRAC fluxes follow a power-law, 
which is consistent with a very dusty active galactic 
nucleus (AGN). 

AzTEC J100006. 5+023837.7 (AzTEC6) - AzTEC6 
is detected at 6.5a by the SMA. It has no candidate 
optical, radio, or 24 /im counterparts, but is detected by 
IRAC at 3.6 and 4.5/im. 

AzTEC 100018.1+024830.5 (AzTEC7)- AzTEC7 is 
detected at 8a by the SMA. Like AzTEC5, there is a 
radio counterpart (i^20cm = 196 ± 61 /iJy), and it is 
detected by IRAC at 3.6, 4.5, 5.8, and 8.0 /im , and 
by MIPS at 24/im. Its SED peaks at 5.8/im, and is 
very bright at 24/.tm, which is consistent with a z ^ 2.5 
starburst. There is also an optical counterpart in the 
ACS imaging with a disturbed morphology, reminiscent 
of a merging system. 
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TABLE 2 

Photometry of SMA/AzTEC Sources 





FuOOtim 

(mJy) 


(mJy) 


rpa 
3.6/j,m 

(MJy) 


(MJy) 


-pa 
b.8fj,m 

(MJy) 


(MJy) 


pb 

(MJy) 


(MJy) 


AzTECl 


10.7 ± 1.3 


15.6 ± 1.1 


4.6 ±1.0 


4.6 ± 1.4 


< 11.2 


17.6 ±8.1 


< 71 


48 ± 14 


AzTEC2'' 


9.0 ±1.3 


12.4 ± 1.0 












52 ± 14 


AzTEC3 


7.6 ±1.2 


8.7 ± 1.5 


3.9 ±1.0 


3.6 ± 1.4 


< 11.2 


< 13.4 


< 71 


< 41 


AzTEC4 


6.8 ±1.3 


14.4 ± 1.9 


4.8 ±1.0 


5.1 ± 1.4 


12.4 ±6.7 


19.6 ±8.1 


< 71 


< 41 


AzTECS 


7.6 ± 1.3 


9.3 ± 1.3 


8.8 ±1.0 


9.0± 1.4 


11.6 ±6.7 


32.7 ±8.1 


164 ± 20 


161 ±35 


AzTECe 


7.9 ± 1.2 


8.6 ± 1.3 


2.4 ± 1.0 


2.4± 1.4 


< 11.2 


< 13.4 


< 71 


< 41 


AzTECT 


8.3 ± 1.4 


12.0 ± 1.5 


52.1 ± 1.0 


52.1 ± 1.4 


80.6 ±6.7 


63.4 ±8.1 


550 ± 20 


196 ± 61 



^ Fluxes are measured in a 3" aperture. Errors and flux limits are the 3cr rms and 5cr rms fluctuation within 
that aperture respectively. Aperture corrections were done to the IRAC calibration radius of 12". 
^ Fluxes were measured in a 5" aperture. Errors are the 3cr rms fluctuations in that aperture, and flux limits 
are the 5a sensitivity from [Sanders et al.l 1 120071 ). Aperture corrections were done to the MIPS calibration 
radius of 35". 

-p20 em measurements have been corrected by a factor of 1.15-1.2 for bandwidth smearing fsee lBondi et all 
I2007IV Flux hmits are at 3ct. 

IRAC and MIPS are confused with a bright foreground object. 



5. DISCUSSION 

Astrometry with the SMA highhghts the unique power 
of this instrument; secure multi- wavelength counterparts 
for many of the targets could only be identified via in- 
terferometric imaging. We find that, while there is al- 
ways - with the exception of the highly confused case 
of AzTEC2 - an IRAC counterpart coincident with the 
SMA position, there are often several 24 /im sources 
within the AzTEC beam (see Figure[T]) that are not asso- 
ciated with the submillimeter emission. This is contrary 
to the prevailing wisdom, in which radio-dim SMGs, like 
their radio-bright counterparts, are associated with red- 
shifted strong polycyclic aromatic hydrocarbon (PAH) 
emission features in the 24 fj.m band. For all five of the 
radio-dim sources (AzTECl-4 and AzTEC6) there are 
proximate 24 /im sources that, if selected, would lead to 
misidentification of multiwavelength counterparts to the 
submillimeter source (see Figure [ij. AzTEC5 has two 
potential radio counterparts within the AzTEC beam, 
both of which have strong 24 fj,m emission. In this case 
as well, only the SMA can unambiguously identify the 
correct counterpart. 

Only now, with these counterparts, can we look for 
clues to the nature of the radio-dim SMC popula- 
tion. It has been suggested that the ratio of the 
submillimeter to radio flux (-F8.';n um/-F2nnm) is a po- 
tentially useful redshift indicator (iCarilli fc YunI Il999l : 
lYun fc Cariliil [200l lAretxaea et al.ll2007D . In"Figure [H 
we plot this ratio for both the SMA/ AzTEC sources 
and for radio-b right SMGs with opti cal spectroscopic 
redshifts (COS: IChapman et all I2005D . The location 
of the radio-dim sources above the locus of points 
from COS is, as pre v ious authors have speculated (e.g.. 
Carilli fc YunI [1991 iChaoman et all [2005t iPope et all 
20061 : llvison et all l2007ir consistent with either colder 
dust temperatures or a higher average/median redshift 
than the COS sample. 

The IRAC and MIPS counterparts for the radio-dim 
SMA/ AzTEC sources support this hypothesis. Given the 
1.1 mm selection function and cosmic volume probed, 
we would expect a sample with Fnoo/jm > 6.5 mJy to 
be dominated by objects at z > 1 with total infrared 
luminosities on the order of 1O^^^^^L0. For such sys- 



tems, the observed 3.6 /im flux should be a strong func- 
tion of redshift. In Figure [3l we show 3.6 /im fluxes 
for the SMA/ AzTEC sources, compared to COS sources 
in the Hubble Deep Field North (HDFN) and SSA22. 
The radio-dim sources have systematically lower fluxes 
which, assuming an Arp 220 model, suggests that they 
may lie at higher redshift. Furthermore, these same 
radio-dim sources do not have 24 /im detections, indicat- 
ing that either the 7.7 /xm PAH emission features have 
been redshifted out of, or the 9.7 /im silicate absorption 
feature is in the 24 /im MIPS band, which suggest ei- 
ther z > 3 or z ^ l.S respectively. However, as a caveat 
wc note that it is not impossible that they are intrinsi- 
cally different objects than Arp 220, with fainter PAHs 
or different submillimeter-to-radio ratios. At the same 
time, radio-bright SMA/ AzTEC sources are consistent 
with the population observed by COS. 

Therefore, the infrared and radio properties of radio- 
dim SMA/ AzTEC sources in a 1.1 mm flux-limited sam- 
ple suggest that they lie at high redshift z ^ 3. Further- 
more, their ubiquity in both this (70% of sources with 
S/N > S at 1.1 m m) and other (S0% of sources with 
S/N > 4 at 1.2 mm: lBertoldi et al.ll200% millimeter sur- 
veys further suggest that they contribute significantly 
to the observed millimeter/submillimeter number counts, 
and that the median redshift of z ~ 2.S inferred by COS 
using assumed radio counterparts is a lower limit. A de- 
tailed discussion of the properties of these systems awaits 
a full SED analysis, which we postpone to a future pa- 
per. However, assuming an Arp 220 model, the observed 
submillimeter flux of these sources implies very high to- 
tal infrared luminosities of L(8 — 1000/im) <; 10^^ Lq - 
or ^ S X 10^^ Lq assuming a Mrk 231 model - which is 
similar to h igh-redshift hyperl uminous infrared galaxies 
(HvLIR.Gs: [Huang et al.ll200l . The presence of a sig- 
nificant population of these objects has important con- 
sequences for models of hierarchical galaxy formation, 
which are only beginning to account for such systems at 
later epochs (z ~ 2; Baugh et al. 200S). It is also curious 
that a majority (60 or 70%) of the most luminous sources 
in the Az TEC/COSMOS ca talog are radio-dim; as first 
noted bv llvison et al.l (j2002f ). the most luminous SMGs 
may have a higher average/median redshift. 
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Fig. 1. — Stamp images for the SMA/AzTEC sources (top to bottom AzTECl-7) for (left to right) SMA (890/xm), VLA radio continuum 
(20cm), MIPS Channel 1 (24/^m), IRAC Channel 1 (3.6M m). and ACS (i-b and; O-S^tm) imaging data. Overlayed in red on the SMA 
image are contours at 3a, 4(t,... from AzTEC imaging data l lScott et al.ll200"^ . The red circles in the remaining stamps have a radius 2", 
corresponding to twice the FWHM of the SMA beam, at the SMA position. Two sources (AzTEC 1 & 7) have secure optical counterparts 
in the ACS images, while AzTEC 3 & 6 have candidate optical counterparts that are a potential foreground object (see §|4]l and outside 
the SMA beam respectively. Each stamp image is 37" X 27", with the exception of the ACS stamps which are 15" X 11". 

our objects (Fsgo/jm ~ 10 mJy) implies dust masses of 
order 0.6 — 3 x 10^ Mp. at the observed time. If dust 



Furthermore, such a population of massive, dusty star- 
bursts at z ^ 3 constrains models of dust production, 
given the limited look-back ti me since the formation of 
the first stars at z « 20 - 30 (jBromm &: LarsonI |2004 
and references therein). The dust mass corresponding 
to the observed thermal emission is approximately equal 
to Aid ~ L^/4:TrK^Bi,{Td), where L^, is the observed lu- 
minosity at a given rest-frame frequency z/, is the 
dust opacity at that frequency, and B^{Td) is the black- 
body emission at the effective du st temperature Td- For 
redshifts of z > 3, assuming the IWeingartner fc Draind 
(j2001i) Milky Way dust opacity, dust temperature Td = 
45 - 70 K, and a flat ACDM cosmology, we find that 
the observed 345 GHz (rest-frame > 1380 GHz) flux of 



10 at the observed time, 
production is dominated by evolved, post- main sequence 
stars with ages ^ 1 Gyr, as it is locally (|Gehrj 119891 : 
iMarchenkol |2006( ) , this requires a dust production rate 
of Md ^ 0.7 — 3.4 Af0 yr~^ over the same dust tem- 
perature range. Or, if supcrnovae are also significant 
contributors to du st production at high rcdshift (e.g., 
iDunne et ani2003fl. as may b e the case in high redshift 
quasars (jMaiolino et al.l[2004[ ). then the required rate of 
production is lower by ~ 50%. 

Finally, SMA imaging in combination with a red- 
shift constraint allows us to place limits on the spa- 
tial extent of the submillimeter continuum. All seven 
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Fig. 2. — The change with redshift of the ratio of the submil- 
hm eter (850Atm) versus radio (20cm) continuum emission in SMGs 
(see lCarilli &: Yunll999l). Radio-b right SMGs with optical spectro- 
scopic redshifts from lChapman "et al. (2005) arc shown as open tri- 
angles, as compared to radio— dim (dashed lines) and radio-bright 
(dash-dot lines) SMA/AzTEC sources (see §|4]for abbreviations). 
The filled triangle represents a rough redshift of z > 4 for AzTECl 
from the B-band dropout nature of its optical counterpart. The 
solid line is a model track for Arp 220. SMA flux measurements at 
890 fim were corrected to 850 ^m using the -Fggo/jm /^i 100^771 ratio 
from SMA and AzTEC, a < 15% correction for all the sources. 




Redshift 

Fig. 3. — IRAC 3.6 fim fluxes versus re dshift for radio— brigh t 
SMGs with optical redshifts (open triangles: IChapman et al.l2005l) . 
as compared to radio-dim (dashed lines) and radio-bright (dash- 
dot lines) SMA/AzTEC sources (see §|3]for abbreviations). The 
filled triangle represents a rough redshift of 2 > 4 for AzTECl 
from the B-band dropout nature of its optical counterpart. For 
comparison, we include model tracks for Arp 220 with total lumi- 
nosities of lO^iL© (a LIRG), lO^^L© (a ULIRG), and IO^'-^Lq (a 
HyLIRG). 

of the sources are compact single sources; the real visi- 
bility amplitudes indicate that they are unresolved out 
to the longest baselines, from which we infer a maxi- 
mum angular size of ^ 1.2". This is particularly in- 
teresting for this sample of the brightest sources in the 
AzTEC/COSMOS catalog because it rules out blends of 
multiple fainter sources as a significant contributor to the 
upper end of the observed SMG luminosity function. It 
also agrees with previous interferometric measurements 
of th e angular exte n t of t h e millimeter ( Downes et al. 
19991: iGen zel ct al. 2003; Downe s fc SolomonI 120031 : 
Kneib et al..,2005. : .Tacconi et al..,2006fr and submillimeter 



(|Iono et al.ll2006f) emission from SMGs, and marginally 
with those of the radio continuum ([Chapman et al.) 
l2004f) . Assuming a flat ACDM cosmology, these 
angular constraints correspond to a physical scale 
for the submillimeter continuum of 10 h^^ kpc at 
z ~ 2 and 8 kpc at z ^4. These size 
scales are consistent with far-infrared continuum emis- 
sion associated with a m erger driven starburst (e.g., 
iMihos fc HernquistI Il994( l analogou s to local luminous 
and u ltralumin ous infrared ga l axies dP^ownes fc SolomonI 
119981: [Sakamoto et all Il999l 120061: llono eTal. 200 7,), 
and are potentiall y in conflict with cool extended cir- 
rus dust models (Efstathiou fc Rowan- RobinsonI l2003t 



iKaviani et al.|[2003 ) and a monolithic collapse scenario. 



6. CONCLUSION 

We use the SMA to follow-up the bright est millime- 
ter sources in the AzTEC/COSMOS survey (jScott et all 
I200I . All seven sources, including five radio-dim 
{F20cm < 60 ^Jy) SMGs, are detected at high signifi- 
cance (> 6<j) with derived positions accurate to ^ 0.2". 
All seven of the sources, with the possible exception of 
one highly confused case, are coincident with IRAC de- 
tections, and all but two are optical dropouts. We find 
that the radio-dim SMGs in our sample have system- 
atically higher submillimeter-to-radio ratios and lower 
IRAC 3.6-8.0 /ini fluxes than radio-bright sources, and 
are not detected at 24/im. This, in combination with size 
constrains from the imaging data, suggests that radio- 
dim SMGs represent a population of very dusty z ^ 3 
starbursts with physical scales similar to local ULIRGs. 
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